Unraveling of free carrier absorption for terahertz radiation in heterostructures 
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The relation between free carrier absorption and intersubband transitions in semiconductor het- 
erostructures is resolved by comparing a sequence of structures. Our numerical and analytical 
results show how free carrier absorption evolves from the intersubband transitions in the limit of an 
infinite number of wells with vanishing barrier width. It is explicitly shown that the integral of the 
absorption over frequency matches the value obtained by the f-sum rule. This shows that a proper 
treatment of intersubband transitions is fully sufficient to simulate the entire electronic absorption 
in heterostructure THz devices. 



I. INTRODUCTION 



are used). The matrix element 



The absorption of electromagnetic radiation due to the 
interaction with electrons in bulk crystals is essentially 
determined by two distinct effects: (i) The free carrier 
absorption (FCA), which is directly related to the elec- 
trical conductivity and drops with frequency on the scale 
of the inverse scattering time, (ii) Interband transitions, 
which are typically described via the dipole moments in- 
duced by the coupling between states in different bands. 
For most crystals these transition energies are of the or- 
der of eV and thus this dominates the response around 
the optical spectrum. In addition to these electronic fea- 
tures, optical phonons provide absorption in the far in- 
frared region, which is not addressed here. 

Semiconductor heterostructures provide an additional 
effective potential for the electron in the conduction band 
causing a further quantization of the electronic states in 
the growth direction (denoted by z). Taking into account 
the degrees of freedom for motion in the x, y plane, this 
establishes subbands within the conduction band. Com- 
monly, the absorption between these subbands is treated 
analogously to the interband transitions in bulk crystals. 
The standard treatment relies on the envelope functions 
(Pv(z) for the subbands v with energies E v and areal elec- 
tron densities n u using expressions for the absorption co- 
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where < E v and the counter-rotating terms are ne- 
glected. Here e is the elementary charge, y/e is the refrac- 
tive index, and eo is the vacuum permeability (SI units 
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describes the coupling strength. Throughout this work 
we assume the polarization of the electric field to point 
in z-direction and that the wave propagates in a waveg- 
uide of effective thickness L z which is filled by the (lay- 
ered) semiconductor material. This scheme is also rou- 
tinely applied for the calculation of the gain spectrum 
of quantum cascade lasers (QCLs). 3 In this context the 
broadening V can be cither added in a phcnomcnological 
way 4 or by detailed calculations, see, e.g., Ref. 5. It can 
also be seen as a limiting case of a full quantum kinetic 
calculation. 6 

While the conventional treatment of intersubband 
transitions is well accepted for transitions in the infrared, 
this approach is less obvious for THz systems, which have 
become of high interest. 7,8 Here, FCA-related features 
might turn up as a strong competing mechanism to the 
intersubband gain transition in analogy to the bulk case 
where both FCA and interband transitions occur as sep- 
arate processes. In order to demonstrate the potential 
relevance, we consider the standard expression for FCA 
in bulk systems 9 
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where m c is the effective mass, n c the volume density 
of electrons in the conduction band, and r is the scat- 
tering time. As an example for GaAs with a doping of 
1 x 10 16 /cm 3 and r = 0.2 ps (corresponding to a mobility 
of 6000 cm 2 /Vs at 300K 10 ) one obtains a = 120/cm for a 
frequency lj/2t: = 2 THz. This is larger than typical gain 
coefficients in THz quantum cascade lasers. 11-13 Thus, 
bulk FCA would provide a strong obstacle in achieving 
lasing in such structures and its proper treatment in het- 
erostructures is of crucial importance for the description 
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TABLE I. Key parameters obtained for the different SLs. The 
effective mass vhsl is taken for the lowest miniband at k — 
in the SL direction. 



barrier width (nm) 


0.5 


1.5 


2.5 


3.5 


miniband width (meV) 


42.7 


25.4 


15.6 


10.1 


effective mass msh/Trie 


0.071 


0.090 


0.125 


0.178 


a at 300 K (A/Vcm) 


11.2 


6.4 


2.4 


0.8 


g at 77 K (A/Vcm) 


17 


14.9 


7.3 


2.9 



of QCLs or other THz hctcrostructurc devices. (For a 
typical infrared laser, in contrast, it was shown that FCA 
in the cascade structure does not play a role. 14 ) In Ref. 4 
FCA was only considered in the waveguide layers but 
not the QCL structure itself, where the absorption was 
determined by intersubband transitions. Furthermore, 
in Ref. 15 it was shown that processes as described by 
Eq. (1) dominate the absorption of light (with z-polarized 
electric field) for quantum wells. 

In this context the question arises how such a treat- 
ment based on intersubband transitions is related to the 
FCA in the bulk. Is FCA related to the seemingly dom- 
inating intersubband processes or does it stem from fur- 
ther processes not identified yet? In the latter case, such 
processes could strongly alter the THz performance of 
heterostructure devices. In order to shed light on this im- 
portant issue we present a detailed study on the unfolding 
of FCA starting from different types of heterostructure. 
Our main conclusion is that the absorption due to inter- 
subband transitions evolves into the bulk FCA for vanish- 
ing barrier widths. This shows that a proper treatment 
of intersubband transitions provides a complete descrip- 
tion of gain and absorption processes in heterostructure 
devices. 



II. FROM SUPERLATTICE TO BULK 

We consider four GaAs-Alo.3Gao.7As superlattices 16 
(SLs) with constant period d = 10 nm. The barrier 
width is set equal to 0.5 nm, 1.5 nm, 2.5 nm, and 
3.5 nm, respectively, and a homogeneous doping with 
n c = 6 x 10 16 /cm 3 is used. The sample with the 2.5 nm 
barrier has been investigated in Ref. 17, which motivates 
our choice. Fig. 1(a) shows the calculated minibands 
assuming effective masses of 0.067m e and 0.0919m e for 
GaAs and Alo.3Gao.7As, respectively, where m e is the 
free electron mass, as well as a conduction band offset 
of 276 meV. 18 Further information on the structures is 
given in table I. 

Here the zero-field conductivity do is evaluated from 
the noncquilibrium Green's function (NEGF) model fol- 
lowing Ref. 19, which includes scattering processes from 
phonons, impurities, interface roughness (with an aver- 
age height of one monolayer and a length correlation 
of 10 nm), and alloy disorder in an approximate way. 
This program also calculates the absorption in linear re- 
sponse to the optical field 6 as given in Fig. 1(b). Using 
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FIG. 1. (Color online) (a) Lowest three minibands for the 
SLs together with the dispersion of bulk GaAs (dotted line) 
neglecting nonparabolicity. (b) Absorption at 300 K for 
the SLs, calculated by the NEGF model together with the 
Drude expression (3) for bulk absorption (dotted line) us- 
ing m c = ttisl = 0.071m e and t — 46 fs. The inset shows 
the NEGF calculations (symbols) at 77 K together with the 
corresponding result of Eq. (5) using ao from table I and 
T m = 70/100/110/100 fs for the sample with the barrier width 
of 0.5/1.5/2.5/3.5 nm, respectively (lines). 



do ~ n c e 2 r/mgL the conductivity for the 0.5 nm bar- 
rier structure provides a scattering time of t = 46 fs. 
This value agrees roughly with the momentum scatter- 
ing rate l/r m = 29/ps (which is the sum of the elas- 
tic and inelastic scattering rate 20 ) extracted from several 
highly doped GaAs/AlAs SLs with narrow barriers at 
room temperature. 21 This value is much smaller than the 
bulk scattering time of 0.2 ps, as scattering is enhanced 
due to the presence of rough interfaces in all SLs (which 
are particular strong scatterers for small barrier widths, 
when the wave functions highly penetrate through the 
barriers). In addition, the assumption of a constant scat- 
tering time is only expected to be of semi-quantitative 
nature, the same holds for the approximations in matrix 
elements used. (For a more detailed treatment of rough- 
ness scattering in thin barriers, see Ref. 22.) Using r = 46 
fs, the Drude expression (3) fits the absorption quite well, 
demonstrating, that these small barriers actually provide 
almost the bulk free carrier absorption behavior. 

With increasing barrier thickness the conductivity be- 
comes smaller due to the reduced coupling between the 
quantum wells. Accordingly, there is a decrease in the 
low frequency absorption 
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as it follows from electrodynamics. 23 Here our numerical 
calculations are in full agreement, as we do not employ 
the rotating wave approximation and include broaden- 
ing in a fully consistent way. Furthermore, for thicker 
barriers, the absorption between the minibands becomes 
more prominent and thus the absorption increases close 
to the photon energy required to overcome the gap be- 
tween the first and the second miniband, as indicated by 
the arrows in Fig. 1(b). The shift of the peak positions 
with respect to the minigaps can be related to scattering 
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induced level shifts. For the 2.5 nm barrier the results 
are in good agreement with the measurements reported 
in Ref. 17. The onset of absorption around 100 meV 
is slightly sharper in the experiment, which may be at- 
tributed to less rough interfaces or to the limited accu- 
racy of the various approximations used for the scattering 
potentials. 

For SLs the absorption can be understood within the 
common miniband picture. For low frequencies intra- 
miniband processes dominate, which are easily treated in 
semiclassical transport models providing for zero electric 
field 24 ' 25 
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This behavior was experimentally observed in Refs. 26 



and 27. Here, ctq 
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widths. With decreasing miniband width, the increase 
of msL reduces errj- An even stronger reduction arises, 
if the miniband width drops below either ksT or the 
Fermi energy, sec Ref. 20 for details. 28 For all supcrlat- 
tice structures studied by our NEGF model, we found 
good agreement with (5) for low frequencies. Some ex- 
amples are shown in the inset of Fig. 1(b). As a further 
example, the calculated absorption spectrum at 65 K for 
the structure of Ref. 27 can be fitted by r m = 0.16ps (not 
shown here). This is in good agreement with the experi- 
mental value of 0.18 ps, which demonstrates the quality 
of the NEGF approach. 

For higher frequencies, transitions between the mini- 
bands can describe the absorption between 60 and 200 
meV very well. See, e.g., the results of the calculations 
in Ref. 17, which fully agree with our more sophisticated 
NEGF approach. 

We conclude that the absorption of SLs at zero bias 
can be well described by the Drudc-like miniband conduc- 
tion result (5) for low frequencies and by common inter- 
miniband transitions for higher frequencies. As shown in 
Fig. 1(b), the combination of both features evolves into 
the bulk FCA (3) if the barrier width becomes small. 



III. FROM MULTIPLE WELL TO 
SUPERLATTICE 

Now we want to study, how the SL absorption arises 
from the behavior of systems containing few wells, which 
show distinct absorption peaks between discrete levels. 
Fig. 2 shows the absorption for multi-quantum-wcll struc- 
tures, as presented in Fig. 2(a) for the case of two wells. 
Here, all parameters correspond to the SL with a 1.5 nm 
barrier discussed above. For the double well structure, 
essentially the two lowest subbands are occupied in ther- 
mal equilibrium, and one observes clear absorption peaks 
corresponding to the separations between the subbands, 
see Fig. 2(b). As the dipole matrix element (2) vanishes 
for equal parity of the states, not all possible transitions 
are visible. The observed peak structure can be directly 
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FIG. 2. (Color online) (a) A double quantum well (well widths 
8.5 nm, barrier width 1.5 nm) with its lowest eigenstates. 
Dashed and dot-dashed lines refer to symmetric and antisym- 
metric states, respectively. The arrows depict the transitions 
associated with the peaks in the absorption spectrum, (b) 
Absorption spectrum calculated by the NEGF model for the 
double quantum well and the corresponding SL. (c) Evolu- 
tion of the absorption for 2, 3, and 4 wells with the same 
parameters as the double well from (a). In order to obtain 
absorption in the entire waveguide, it is assumed that the 
multi-quantum-well structure is periodically repeated with a 
separation by a 7.5 nm barrier. All calculations are done at 
T = 300 K. 



described by the standard intcrsubband expressions (1). 
Furthermore, there is zero absorption in the limit of zero 
frequency as no dc current along the structure is possible, 
compare Eq. (4). 

With increasing well numbers, the peaks III and IV of 
the double well split up and form the continuous absorp- 
tion between 60 and 200 meV due to the transitions be- 
tween the first and the second SL miniband, see Fig. 2(c). 
While this is quite expected, peak I does not show any 
clear splitting, but shifts to lower frequencies, approach- 
ing the intra-miniband absorption. This behavior can be 
understood by a detailed study of the multi-quantum- 
well eigenstates. Here, a tight binding model for N wells 
with next neighbor coupling T± shows (see the appendix 
A for details): (i) There are N eigenstates, labeled by an 
index v according to their energy E v . Here E u+ \ — E v 
is of the order of 4|Ti|/7V. (ii) The matrix element 
from Eq. (2) is small unless for neighboring states, i.e. 
fx = v ± 1. Thus, the transitions between neighboring 
states dominate, explaining the strong absorption around 
hw « 4]Ti]/JV visible in Fig. 2(c), where 4|Ti| essentially 
corresponds to the miniband width of the infinite struc- 
ture. Together with a tail at higher frequencies due to 
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broadening of these transitions this explains the appear- 
ance of the Drude-like miniband absorption for the SL 
in the limit of large N. For oj = the evolution is not 
smooth as any finite sequence of quantum wells has a zero 
dc conductivity in contrast to an infinite SL and thus the 
absorption must vanish according to Eq. (4). 



IV. THE INTEGRATED ABSORPTION 

Summing over all possible intersubband transi- 
tions (1), we obtain the total absorption ais(oj) = 
Smj/ Q^vi^QiEv — Efj). Here the discrete index v runs 
over all (infinitely many) eigenstates of the heterostruc- 
ture of finite length, including states which correspond 
to unbounded states with energies far above the barrier 
potential. Integrating over all frequencies provides 
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under the assumption E v — E^ ^> T - otherwise the 
counter-rotating terms become of relevance, which had 
been neglected here. In appendix B we show that the 
same integral relation is more generally obtained for ar- 
bitrary level spacings E v — E^ within our NEGF model, 
which also covers dispersive gain. 29 ' 30 

Following Ref. 31, Eq. (6) can be simplified by 
the Thomas-Reiche-Kuhn sum rule 32 (also called f-sum 
rule 33 ) which reads for a parabolic band with effective 
mass m c 
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and provides the integrated absorption 
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where rt av = ^2 n^/L z is the average three-dimensional 
carrier density in the waveguide. 

For a bulk semiconductor, the free carrier absorption 
(3) provides after integration over energy 
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which fully agrees with the intersubband result (7) for 
equal total densities n c = n av . Thus the total FCA 
in a bulk semiconductor equals the total intersubband 
absorption within the conduction band for a finite het- 
erostructure of finite length, which shows the direct rela- 
tion between these. More generally, Eqs. (7,8) establish 



a general rule for the integrated absorption within the 
conduction band of a semiconductor under conditions, 
where the approximation of a constant effective mass is 
justified. In this context superlattices appear as an inter- 
mediate case, where the inter-miniband absorption and 
the Drude-like intra-miniband absorption add up to the 
full result. 31 

Our numerical data in Fig. 1(b) exhibit the inte- 
grated absorption ft J du a(oj) = 25.7(±0.3)eV/cm for 
all curves. The data from Fig. 2 provide ft J dtu a(oS) = 
N ^ Q 6 25.6(±0.2) eV/cm, where the additional factor 
takes into account the undoped region of 6 nm between 
adjacent multiple quantum wells (N is the number of 
well/barrier combinations with a length of 10 nm each). 
These values are slightly below the value of 27.3 eV/cm 
given by Eq. (8) using the GaAs effective mass. This mi- 
nor discrepancy of less than 7% can be easily attributed 
to some absorption at higher frequencies and the impact 
of the barrier material with a larger mass. 34 We conclude, 
that the absorption obtained by our NEGF code is in ex- 
cellent agreement with the rule (7). 

More generally the effect of the semiconductor het- 
erostructure can be understood as shifting the absorp- 
tion strength within the frequency space, as explicitly 
demonstrated by our calculations. This perception has 
actually been used in the design of QCL structures, where 
the unavoidable free carrier absorption is deflected from 
the frequency region of operation by a proper choice of 
heterostructures 35 , see, e.g., Ref. 36. 



V. CONCLUSION 

We demonstrated how the common bulk free carrier 
absorption evolves from standard intersubband absorp- 
tion in heterostructures for electromagnetic waves with 
an electric field pointing in growth direction. Here the 
well-studied SL absorption constitutes an intermediate 
case, which can be entirely understood on the basis of 
common intersubband absorption processes in the limit 
of a growing number of quantum wells. For decreasing 
SL barrier width the combination of inter- and intra- 
miniband absorption evolves into the standard FCA of 
the bulk crystal. This behavior reflects a redistribution 
of absorption strength, while the integrated absorption is 
constant. The most relevant consequence is that there is 
no need to bother about any additional FCA-related ab- 
sorption processes, provided all intersubband transitions 
are properly taken into account. A consistency check 
for the calculated gain/absorption spectrum is whether 
Eq. (4) is satisfied in the low frequency limit and the 
integrated absorption matches Eqs. (7,8). 
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Appendix A: Analytical calculation for coupled wells 



We consider a multi-quantum-well structure with N 
wells centered at z = nd, where n = 1,2, ...N. The 
ground state of the isolated well n has the wavcfunction 
\£" g (z — nd) and the energy E g . Restricting to a nearest 
neighbor coupling T\ (which is negative for the lowest 
subband), the eigenenergies are 



E,, — E„ 



2T X cos 



AT + 1 



fori/ = 1,2,... JV (Al) 



and the eigenstates read <-p v {z) = a^^ g (z — nd) with 
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If the overlap between the states in different wells is neg- 
ligible, i.e. J dz^ g (z — n' 'd) g (z — nd) « nd5 nn >, we find 



z^ v = ~^2 n nddn ctn j which can be directly evaluated. If 
z/ — is even we find = S fJiM (N + l)d/2 as both states 
have the same parity with respect to z = (N + l)d/2. 



For odd v — ii, some algebra yields 
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For fi/v we thus have: 
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z^jy = O | ^ — _^ 1 for odd and large [y — fx) 



As the square of z^ v enters the absorption (1), it becomes 
clear that the transitions with v = u ± 1 highly domi- 
nate the absorption spectrum. The energy difference of 
the corresponding states (Al) for these transitions is less 
than 2|Ti|7r/(7V + 1) with an average of approximately 
4|Ti|/JV. 



Appendix B: Total absorption with the Green's 
function model 

Here we refer to the formulation of our NEGF model 
as outlined in Rcf. 6. Here gain is evaluated within lin- 
ear response around the stationary state characterized 
by the Green's functions G MI ,(k, E). In order to simplify 
the analysis, nondiagonal G„„(k, 2?) are neglected here - 
they are, however, fully included in our numerical imple- 
mentation. Then the absorption resulting from the pair 
of states fx, v can be written as 



^M = e2( ^f^ |2 |E / ^{G- t (k,E + MG< M (k,iJ) + (5<(k,iJ + MG^(k,s) 

k (Bl) 

- G^(k, E + MG< (k, E) - G< M (k, E + MG^ v (k, E)} 

which is essentially the last equation of the appendix in Ref. 6 with the counter-rotating term added. Inserting the 
spectral function 37 A v (k,E) = T23{G^ /adv (k, E)} and its occupied part 38 A° cc (k,£) = -iG< v (k,E), which is 
assumed to be real, we find 

a ^H- e2( ^f^' 2 |E / ^[^ cc ( k '^( k ^ + M-^ cc (k,^(k,£;-M ^ 

+ A° cc (k, E)A v {k, E - hw) - A° cc (k, E)A^{k, E + hio) 

The terms A° cc (k, E)A V (k, E+hw) — A° cc (k, E)A^(k, E—fno) provide the physical origin of dispersive gain as sketched 
in Refs. 6 and 30. The signs of the counter-rotating terms ^4° cc (k, E)A l/ (k, E - hio) - A° cc (k, E)A^(k, E + Huj) seem 
to contradict our intuition, as the first one appears to relate to emission and the second to absorption. However, in 
this formulation the sign is defined via the difference in energy between the initial and the final state, where only one 
a specific combination is used in the prefactor (E v — E^). 
Using the general relations 

^ dcj[A^k,E + hu) + A l ,(k,E-fiw)} = -j dE'A fl (k,E')=2n/h and - ^ j _A° cc (k, £) = n p 



6 



integration of the terms from Eq. (B2) over frequency provides 

' oW^H = ^ 2 M{E V -E ){n^n v ) 

so that the sum over all different pairs {n,v) equals the second line of Eq. (6). Thus the integrated absorption (7) 
also holds for the more involved absorption terms (B2) of the NEGF model which include the dispersive gain. 

I 
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